56 research outputs found
Three-level BDDC in three dimensions
This is the published version, also available here: http://dx.doi.org/10.1137/050629902.Balancing domain decomposition by constraints (BDDC) methods are nonoverlapping iterative substructuring domain decomposition methods for the solution of large sparse linear algebraic systems arising from the discretization of elliptic boundary value problems. Their coarse problems are given in terms of a small number of continuity constraints for each subdomain, which are enforced across the interface. The coarse problem matrix is generated and factored by a direct solver at the beginning of the computation and it can ultimately become a bottleneck if the number of subdomains is very large. In this paper, two three-level BDDC methods are introduced for solving the coarse problem approximately for problems in three dimensions. This is an extension of previous work for the two-dimensional case. Edge constraints are considered in this work since vertex constraints alone, which work well in two dimensions, result in a noncompetitive algorithm in three dimensions. Some new technical tools are then needed in the analysis and this makes the three-dimensional case more complicated. Estimates of the condition numbers are provided for two three-level BDDC methods, and numerical experiments are also discussed
A Nonoverlapping Domain Decomposition Method for Incompressible Stokes Equations with Continuous Pressures
This is the publisher's version, also available electronically from http://epubs.siam.org/doi/abs/10.1137/120861503A nonoverlapping domain decomposition algorithm is proposed to solve the linear system arising from mixed finite element approximation of incompressible Stokes equations. A continuous finite element space for the pressure is used. In the proposed algorithm, Lagrange multipliers are used to enforce continuity of the velocity component across the subdomain boundary. The continuity of the pressure component is enforced in the primal form, i.e., neighboring subdomains share the same pressure degrees of freedom on the subdomain interface and no Lagrange multipliers are needed. After eliminating all velocity variables and the independent subdomain interior parts of the pressures, a symmetric positive semidefinite linear system for the subdomain boundary pressures and the Lagrange multipliers is formed and solved by a preconditioned conjugate gradient method. A lumped preconditioner is studied and the condition number bound of the preconditioned operator is proved to be independent of the number of subdomains for fixed subdomain problem size. Numerical experiments demonstrate the convergence rate of the proposed algorithm
Limitations of polynomial chaos expansions in the Bayesian solution of inverse problems
Polynomial chaos expansions are used to reduce the computational cost in the
Bayesian solutions of inverse problems by creating a surrogate posterior that
can be evaluated inexpensively. We show, by analysis and example, that when the
data contain significant information beyond what is assumed in the prior, the
surrogate posterior can be very different from the posterior, and the resulting
estimates become inaccurate. One can improve the accuracy by adaptively
increasing the order of the polynomial chaos, but the cost may increase too
fast for this to be cost effective compared to Monte Carlo sampling without a
surrogate posterior
Singular function mortar finite element methods
This is the published version, also available here: http://dx.doi.org/10.2478/cmam-2003-0014.We consider the Poisson equation with Dirichlet boundary conditions on a polygonal domain with one reentrant corner. We introduce new nonconforming finite element discretizations based on mortar techniques and singular functions. The main idea introduced in this paper is the replacement of cut-off functions by mortar element techniques on the boundary of the domain. As advantages, the new discretizations do not require costly numerical integrations and have smaller a priori error estimates and condition numbers. Based on such an approach, we prove optimal accuracy error bounds for the discrete solution. Based on such techniques, we also derive new extraction formulas for the stress intensive factor. We establish optimal accuracy for the computed stress intensive factor. Numerical examples are presented to support our theory
A FETI-DP TYPE DOMAIN DECOMPOSITION ALGORITHM FOR THREE-DIMENSIONAL INCOMPRESSIBLE STOKES EQUATIONS
The FETI-DP (dual-primal finite element tearing and interconnecting) algorithms,
proposed by the authors in [SIAM J. Numer. Anal., 51 (2013), pp. 1235β1253] and [Internat. J.
Numer. Methods Engrg., 94 (2013), pp. 128β149] for solving incompressible Stokes equations, are
extended to three-dimensional problems. A new analysis of the condition number bound for using
the Dirichlet preconditioner is given. The algorithm and analysis are valid for mixed finite
elements with both continuous and discontinuous pressures. An advantage of this new analysis is
that the numerous coarse level velocity components, required in the previous analysis to enforce the
divergence-free subdomain boundary velocity conditions, are no longer needed. This greatly reduces
the size of the coarse level problem in the algorithm, especially for three-dimensional problems. The
coarse level velocity space can be chosen as simple as those coarse spaces for solving scalar elliptic
problems corresponding to each velocity component. Both the Dirichlet and lumped preconditioners
are analyzed using the same framework in this new analysis. Their condition number bounds are
proved to be independent of the number of subdomains for fixed subdomain problem size. Numerical
experiments in both two and three dimensions, using mixed finite elements with both continuous
and discontinuous pressures, demonstrate the convergence rate of the algorithms
- β¦